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Abstract 

Various inflationary scenarios can often be distinguished from one another by looking at the 
squeezed limit behavior of correlation functions. Therefore, it is useful to have a framework de¬ 
signed to study this limit in a more systematic and efficient way. We propose using an expansion 
in terms of weakly coupled super-horizon degrees of freedom, which is argued to generically exist in 
a near de Sitter space-time. The modes have a simple factorized form which leads to factorization 
of the squeezed-limit correlation functions with power-law behavior in fciong/^short- This approach 
reproduces the known results in single-, quasi-single-, and multi-field inflationary models. However, 
it is applicable even if, unlike the above examples, the additional degrees of freedom are not weakly 
coupled at sub-horizon scales. Stronger results are derived in two-field (or sufficiently symmetric 
multi-field) inflationary models. We discuss the observability of the non-Gaussian 3-point function 
in the large-scale structure surveys, and argue that the squeezed limit behavior has a higher de¬ 
tectability chance than equilateral behavior when it scales as (Mong/^short)'^ with A < 1 - where 
local non-Gaussianity corresponds to A = 0. 


1 Introduction 

One of the main goals of modern cosmology is to understand the origin of primordial fluctu¬ 
ations. We are fairly conhdent, given the in-phase acoustic oscillations, that the fluctuations 
have been super-horizon for a while. Inflation seems to be the simplest, and most widely 
accepted, explanation for generating the observed slightly red-tilted and almost Gaussian 
spectrum. However, there are many different inflationary scenarios, and there are alterna¬ 
tives to inflation, which we would like, and potentially be able to distinguish using data 
from observations in the foreseeable future. This program requires an understanding of the 
cosmological observables, and a classihcation of different models in terms of their distinctive 
observational signatures. 

Perhaps the most decisive of these signatures whose detection rules out non-inflationary 
alternatives is the H-type polarization of the cosmic microwave background caused by pri¬ 
mordial tensor modes. However, depending on the energy scale of inflation, the amplitude 
of tensor modes may be too small for the H-modes to be observable, and even if observed. 
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there is still much to learn about inflation. In particular, understanding what other degrees 
of freedom besides the inflaton are active during inflation and what is their spectrum of 
masses and couplings can teach us about the underlying ultraviolet physics. 


It is well known through several examples that the information about these degrees of 
freedom is encoded in the squeezed limit of correlation functions in momentum space, when 
one or a partial sum of momenta becomes “soft”, i.e. much smaller in magnitude than 
the others. Single-field models of inflation satisfy a set of identities relating the squeezed 
limit correlation functions to lower order statistics [1, 2]. In particular, denoting the soft 
momentum by q and factoring out the power spectrum of the soft mode P{q), the zeroth 
and first order terms in q are fully fixed by these identities [3, 4]. The model dependent 
self-interactions (for instance, the equilateral type 3-point correlators) start to contribute 
starting from second order. That is, neglecting the tensorial structure. 




nonuniversal 
single—field 


P{g) 


cc 


( 1 . 1 ) 


On the other hand additional light degrees of freedom during inflation leave characteristically 
different imprints in the squeezed limit correlation functions whose scaling in the ratio of 
the small momentum to the large ones q/k is determined by the ratio of the mass to the 
expansion rate during inflation [m/H) [5, 6]. The purpose of this note is to give a more 
unified treatment of this limit. 


We take an approach similar in spirit to the effective held theory (EFT) of single-field 
inflation [7]. In the simplest single-field, slow-roll models of inflation the background ex¬ 
pansion as well as the dynamics of scalar fluctuations are fully described in terms of one 
fundamental scalar field, the inflaton. However, there exist several models where the back¬ 
ground is not driven by a weakly coupled fundamental held, such as DBI and /c-inflation, 
but the fluctuations have a richer phenomenology (e.g. non-Gaussianity). These models are 
generally characterized by having a lower strong-coupling scale compared to the energy scale 
at which de Sitter isometries are broken [8]. The advantage of the EFT of inflation is to pro¬ 
vide a unifying framework to study the perturbations in these models by means of separating 
the discussion of background evolution from the fluctuations and parameterizing the former 
in a set of time-dependent coefficients. This allows a much more efficient and systematic 
treatment of the fluctuations, and in terms of parameters that are directly determined by 
observations. 


Similarly, the simplest inflationary models with more than one degree of freedom are 
usually described by the addition of new fields that are weakly coupled up to an energy scale 
A ^ H. They can be included in the EFT framework as in [6, 9], and their coupling to 
the inflaton, or their subsequent conversion to adiabatic modes modifies the phenomenology 
of inflation. However, distinct features generally appear in the behavior of the correlation 
functions in the squeezed limit. Away from this limit the shapes of correlation functions 
are expected to be degenerate with those resulting from self-interactions in a single-field 
model. Hence, although the EFT of multi-field inflation allows us to calculate the full 
shape of non-Gaussianity, only a part of it is of special interest. We therefore look for an 
effective description of this part by working directly at the level of the correlation functions 
rather than the Lagrangian. We will argue that such description exists, it is independent of 
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details of interactions and is parameterized in terms of quantities that are directly related 
to observations. 

One can also imagine scenarios in which the additional degrees of freedom are strongly 
coupled below and around the Hubble scale, i.e. A ~ The EFT approach can fail in 
this case, while as we will argue, at super-horizon scales one seems to be able to always 
talk about weakly coupled degrees of freedom. This is a consequence of locality and the 
slow-roll assumption, namely that the time-variation of the parameters of the underlying 
theory is much slower than the expansion rate H. Therefore, one advantage of our approach 
is that it proves many of the known results about the squeezed limit behavior of inflationary 
correlation functions hold more universally, beyond the explicit models in which they were 
originally derived. In particular, the squeezed limit behavior is describable in terms of a 
discrete sum of power laws with either real positive or pairs of complex conjugate powers 
whose real parts are multiples of 3/2. Any deviation (such as negative powers or a continuous 
spectrum) seems to indicate a radically different alternative to inflation.^ 

In what follows, we first investigate the evolution of local operators at super-horizon 
scales in more detail, and argue that they have simple scaling behavior as a function of the 
conformal time rj, up to slow-roll corrections. In many cases of interest they become classical 
and their time and space dependence factorizes: 

a{ri,x) ^ a{x)ri'^'^. (1.2) 


We then show that to describe the squeezed limit behavior of the observed correlation func¬ 
tions, all that is needed is the spectrum of {Ag-}, plus a set of coefficients parameterizing 
the degree of coupling between adiabatic fluctuations and each local operator, e.g. 


(C(a?)C(0)), =^o(a?) + d(0)^<,(a;) + dia{0)Q^{x) + didja{0)C^g^{x) + ■■■ 
+ ^^{0)^a^{x) + adia{0)Cg^{x) + ■■■ 


(1.3) 


At this order the ^ functions are hxed up to numerical coefficients which are the theoretical 
data directly related to the observations. This applies to quasi-single held and multi-held 
models but also when ( is directly coupled to a composite operator like The contribution 
of an operator of dimension Aq > 0 to the squeezed limit scales as 


iCq ■ ■ ■ Ck) Q 

Pig) 


oc 


(1.4) 


The special case in which the perturbations of the intermediate operator do not decay sig- 
nihcantly by expansion (Ao — 0) deserves a more detailed discussion. These huctuations 
are called entropic perturbations, which can convert to adiabatic perturbations at a later 
time during or after the end of inhation. The squeezed limit behavior is again simply related 
to a set of mixing and coupling coefficients. In the case where there is only one entropy 
mode, or there are several identical ones, we can determine the degree to which the observed 
adiabatic huctuations are coming from the initial huctuations of the inhaton, and derive a 

^In many of the existing alternatives to inflation the fluctuating fields see an effective de Sitter background 
metric (see [10] and references therein), and hence our conclusions remain unmodified. 
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new consistency condition satisfied by ratios of squeezed limit correlation functions. In the 
more general case, unitarity imposes an inequality condition on these ratios [11, 12], 

Exceptions to (1.2) occur when there are massive particles of m > 3H/2. This case is 
discussed separately since the corresponding super-horizon modes never become classical; 
they are described by two complex conjugate power laws and the notion of particle remains 
applicable. However, there is still a sense in which the correlation functions factorize in the 
squeezed limit which is analogous to that of scattering amplitudes on single particle poles. 
The squeezed limit of the bispectrum now scales as 

P{q) 

but it also oscillates as a function of the ratio and the angle between the short and long 
momenta. This case is particularly interesting since the squeezed limit correlation functions, 
which in real space describe the correlation between short distance ~ 1/k measurements 
performed at large separation ~ 1/g, have quantum mechanical origin in contrast to all 
previous cases. The oscillations arise from interference between two branches of the wave- 
function [13, 14]. 

Finally, we discuss the observational prospects of these squeezed limit behaviors. It 
is known that local non-Gaussianity (a non-Gaussianity type that can be generated by 
entropy modes) has a higher chance of being detected by observation of “non-local biasing” 
in large-scale structure surveys [15]. This is because the squeezed limit behavior of local non- 
Gaussianity is distinct from what can be caused by gravitational dynamics in the formation of 
structures, which scales as in the soft momentum. In particular, by going to larger scales, 
or equivalently sending g —)■ 0 the contribution of local non-Gaussianity rises compared to the 
gravitational one. A milder enhancement can arise in all of the above scenarios. However, 
due to the hnite size of the universe the amount of signal in the long wavelength modes 
is limited. Through an estimate of the signal-to-noise ratio we will argue that only those 
non-Gaussian 3-point functions that have the squeezed limit scaling below 0{q) are better 
detectable in this limit. 


(1.5) 


2 Super-horizon degrees of freedom 

Let us begin by considering an interacting theory of helds of non-zero mass during inflation. 
We denote them by {ca}, and order them with respect to their mass: rria < ma+i- As 
long as the time variation of coupling coefficients in the interaction Hamiltonian "H is slower 
than the Hubble rate (the slow-roll assumption), it is easy to see that interactions become 
irrelevant for super-horizon modes. Hence, one can choose a basis of free fields. 

Note hrst that if one ignores the interactions among the sub-horizon and super-horizon 
modes and focus only on the latter, the terms with spatial derivative are negligible. Thus the 
equations of motion form a system of ODEs in time. Take the equation(s) for the lightest 
held(s). It is of the general form 

+mlai = Fi{{aa},t), ( 2 . 1 ) 


4 



where Fi is at least of second order in fields or their time-derivatives and it can explicitly 
depend on time through time-dependent coupling coefficients. Over-dot denotes dt- To argue 
that Fi becomes negligible outside of the horizon we need to analyze the linear solutions. Up 
to slow-roll suppressed corrections, the super-horizon free-held solutions in a near de Sitter 
space-time of expansion rate H are of the form 


oc 


(g^)A± 

g3/2 ’ 




3 /g ^ 


( 2 . 2 ) 


where rj = —1/aH is the conformal time and the q dependence is hxed by the canonical 
quantization condition. It follows from the solution (2.2) that for masses below 3H/2 the 
mode functions divide into growing and decaying modes. Therefore, after a few Hubble times 
the held huctuations are dominated by a single growing mode and the spatial and temporal 
dependence of the perturbations factorize, giving 




(2.3) 


up to corrections that decay with a negative power of a. It follows that the held has become 
classical: by locality the equal-time real space commutator is proportional 

to a~^S^{x — x') which implies that the momentum space commutator 

W{q, i), ^{q, i)] oc a~^5^{q + q) (2.4) 


decays much faster than the product (j{q,t)&{q',t) oc . Hence the held and its time 

derivative can be simultaneously known with negligible uncertainty. For heavier helds with 
m > 3H/2 this conclusion does not hold. 

Returning to the main problem, the slow-roll assumption about the explicit time depen¬ 
dence of Fi together with rria > 0 then imply that the interactions on the right-hand side of 
(2.1) decay by a negative power of a, and hence become negligible at late times. That is ai 
becomes free. 

We can now proceed to the next held, say <72 with equation of motion given by 

a~^dt{a^&2) + mla2 = F2{{aa}, t). (2.5) 


Let us focus on one example, for instance, F 2 = Xaf (the generalization to arbitrary in¬ 
teraction is straightforward). If then the above argument can be repeated to 

show that F 2 becomes irrelevant at late times. However, in the opposite regime there can 
be interactions made of lighter helds which are not negligible. Here we can use the factor¬ 
ization property of the solutions for light helds to diagonalize the equation (2.5). Note that 
A~^ > 2A~^ can happen only if A~^ <3/4 for which the above-mentioned dominance of the 
growing mode and factorization of super-horizon solution is valid. We can therefore redehne 
the heavier held by adding to it a particular solution of the equation made of the lighter 
helds. In the above example, take 


<^2 — d2 -|- 


2Ai{2Ai - 3)m+ mf^' 


( 2 . 6 ) 
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where Ai = A“^. The field 5'2 now satisfies a free field equation of mass m 2 , plus corrections 
that are negligible outside of the horizon. This procedure can be continued to obtain a 
complete basis of free fields {da}- The tildes will be dropped in the following. 

So far the interactions with modes of the horizon size or shorter have been ignored. Once 
they are included and averaged over two types of corrections arise. First, their back-reaction 
to the long modes renormalizes the mass and coupling coefficients of the theory (or introduces 
new interactions), and secondly, the short modes can stochastically combine to make up a 
long wavelength mode. The first correction does not introduce any qualitative difference in 
the above argument except that the masses must be replaced by the renormalized ones. On 
the other hand, the stochastic contribution of horizon size modes on a super-horizon mode 
q is suppressed because it requires N ~ (aH/q)^ independent Hubble patches to conspire. 
The contribution to the variance of a is therefore expected to be of order {q/aHY and hence 
negligible [16]. 

Since every perturbative theory of massive fields approaches the free theory outside of 
the horizon, we expect any (possibly strongly coupled) theory which can be continuously 
connected to a perturbative theory by dialing its coupling coefficients to also satisfy this 
property. If there are massless fields that are not protected by a symmetry, the interactions 
do not necessarily decay with a. However, keeping these scalar fields massless seems to 
require considerable fine-tuning. We discard this possibility in what follows. Note that 
the metric fluctuations (both tensor and scalar modes) are massless, but their interactions 
contain sufficient number of derivatives and die off outside the horizon. 

It is important to stress that the above argument does not ensure that the super-horizon 
degrees of freedom continue to exist as weakly coupled fields at horizon and sub-horizon 
scales. However some general statements about correlation functions can be made even 
in these cases. For example, one can imagine that the physics at the sub-horizon scales 
is described by a strongly coupled conformal field theory (CFT).^ Or the compactification 
scale of the string theory may happen to be close to H in which case a weakly coupled 4d 
description might not exist. The other interesting example is models with dissipation [18]. In 
all these and similar scenarios it is hard or impossible to calculate the correlation functions. 
Nevertheless, the results of the following sections hold regardless of that. They are based 
just on the existence of free super-horizon degrees of freedom. 


3 Squeezed limit in a classical background field 


The result of the previous section has important consequences for the squeezed limit correla¬ 
tion functions. In order to see that let us consider a long-wavelength perturbation of a light 
field whose space and time dependence factorizes as argued above, a{x, rj) ~ d'{x)g{ri), where 
g{r]) is the growing solution. This long wavelength mode slightly modifies the background 

^One such example was proposed in [17]. There it was assumed that a sector of theory is described both 
at super- and sub-horizon scales by a CFT that is conformally coupled to the metric. We think one can 
consider a more general case where the coupling of the ultraviolet CFT to the metric is not conformal and 
hence it flows to a QFT with a mass gap of order H on curved background. The above arguments would 
then apply in this case. 
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seen by the short wavelength fluctuations. Hence, in the squeezed limit of cosmological cor¬ 
relation functions, the contribution of a is through the response of shorter wavelength modes 
to this background. More specifically, this happens when the short modes start interacting 
with the long mode after its horizon crossing, or if their memory of earlier interactions is 
erased due to adiabatic evolution. In these cases an effective description can be employed. 
First, note that the short scale physics can only depend on local observables made of the 
long-wavelength mode which can be organized in a derivative expansion. For instance, the 
2 -point correlation functions of shorter wavelength metric fluctuations ( on the background 
can be written as 


(C(a?)C(0)), =^o(aj) + d(0)^<,(a;) + dia{0)CdA^) + + '' ' 

+ + ■■■ 


(3.1) 


where .^o is the correlation function in the absence of the particular long wavelength mode. 
This expansion in terms of local operators is not necessarily valid to all orders in the soft 
momentum q. For now, we are mainly interested in the leading contribution, however this 
point will be discussed in more detail in section 3.1.1. 

Equation (3.1) implies that 


{CqCkiA) — + ■ ■ ■ 


(3.2) 


where PqA) cross correlation Pc^A) ~ iCq^-q) and dots represent other contributions 

to the squeezed limit 3-point function. If the mixing r = Pc_{q)PA) ^ 1) then to 

leading order in r we can write 


{CqCkiCk2) a — 


(Cq^-<?(h0)) (CfclCfc2^q(^0)) 

Pa{q,Vo) 


(3.3) 


where rjQ is a time late enough so that all of the modes are super-horizon. The same formula 
(3.1) can also be used to And the squeezed limit of a four-point function when the sum of 
two momenta is much less than individual momenta (this is sometimes called the counter- 
collinear limit). We get (for ki ^ —^2 and ^ —k^) 


(CfclCfc2Cfc3Cfc4) — ^ ^ P^Pim fe)««(fci)U(fc3) + '" (3.4) 


Note that the background long wavelength mode can be a fluctuation in ( (or the tensor 
modes 7 *^). However, unlike the generic case, now there are consistency conditions on how 
the short modes respond to a long-wavelength metric fluctuation. Some linear superposi¬ 
tions of these fluctuations are locally indistinguishable from coordinate transformations and 
the response of the short modes is completely fixed by symmetries [4]. Model dependence 
starts to show up at order didj( which is the reason for the squeezed limit behavior of 
equilateral-type shapes (1.1). In the absence of other super-horizon degrees of freedom the 
above constraints lead to consistency conditions on the full squeezed limit correlation func¬ 
tions, which hold regardless of the details of sub-horizon physics. For instance, an explicit 
derivation when the sub-horizon dynamics is effectively dissipative can be found in [19]. 
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In the presence of other degrees of freedom, two cases are of special interest and will be 
discussed separately in the following sections. First, when a has mass comparable to H (or 
more precisely w?/ 3> 1 — n^, the scalar tilt). In this case a fluctuations gradually decay 
outside the horizon, and the long-short correlation (3.1) forms around the horizon crossing 
time of the short modes during inflation. The second case is when m/H 1. In this case 
(3.1) can form simultaneously for all k, long after they have all crossed the horizon. We 
denote these nearly massless fluctuations as entropy modes. 

3.1 Massive fields with 1 — rtg <C m?/H‘^ <9/4 

Neglecting slow-roll corrections, the growing mode of a held with mass m < 3H/2 at super¬ 
horizon scales behaves as 

(^qiv) - oc (3.5) 

Since the correlation (3.1) forms around the horizon crossing time rj ~ —1/fc, the factor 
translates into k~^'^. On the other hand, the scale invariance of de Sitter (i.e. the symmetry 
under {ri,x) —)■ {\ri,\x)) implies that 


{Cqa-q{rio)y oc —. (3.6) 

From (3.5) we see that Vo) oc {qVo)'^^'" Ip’ up to normalization constants and {aq{rjo)^-q) 
Pa{q, Vo)/Vo'' ■ Hence, the q dependence of (3.2) and (3.3) are fully hxed and using the scale 
invariance we can write 

(C,&.Cb),«/>c(9)fcW(f)"'". (3.7) 

where in this approximation P(^{q) cx q~^. It is seen that a wide range of squeezed limit 
behaviors with 0 < A“ < 3/2 result depending on the mass of a. The quasi-single held 
scenario [20] is the perturbative realization of this possibility. A review of the perturbative 
calculation is provided in appendix A.l. 

Note that there can also be a mixing between C, and a composite operator made of a. 
For instance, if 

fc^-(9) = (C<(9))V0 (3.8) 

then there will be a contribution like (3.7) but with A“ —)■ 2A“. More generally for any 
local operator O of dimension Ao < 2, we have the squeezed limit behavior {q/k)^° that 
decays slower than the equilateral-type contributions. This is explicitly demonstrated in a 
perturbative example in appendix A.2. Note that the whole range between the local and 
equilateral type non-Gaussianity is covered in these models. This type of correlation between 
the metric huctuation ^ and operators constructed out of super-horizon degrees of freedom 
can naturally arise, for instance via a coupling 



3.1.1 Squeezed limit and operator product expansion 

Before discussing the case of entropic fluctuations, let us make some comments which cast 
doubt on any general connection between squeezed limit correlation functions during infla¬ 
tion (or a near de Sitter phase) and operator product expansion (OPE) in a hypothetical 
conformal held theory in one lower dimension (hints that such a connection might exist can 
be found e.g. in [14, 21]). The presence of such a connection is motivated by the fact that 
the action of de Sitter isometries on super-horizon helds can be identihed with the action 
of conformal group on some dual operators in a CFT living on a spatial slice. This can 
best be seen at the level of the wavefunction of the universe. The wavefunction is most 
naturally expressed for helds of dimension lower than 3/2 in terms of sufficiently smooth 
(super-horizon) held conhgurations as \h[{99},77]. This is the analog of coordinate basis for 
particle quantum mechanics but now with one real variable for each iflx), or for 

each 93 (g) with q The expectation values are calculated by integrating over all of 

these helds 

{f’l{vO)}l> = / r'^|4'|{v}..)]pFl{v}|. (3.9) 

Note that 93 on the r.h.s. is an integration variable and the time-dependence is fully encoded 
in T. 

The coefficient of the term 931 ■ ■ ■ 93 ^ in log T have the same structure as an n-point 
correlation function of CFT operators (Oi ■ ■ ■ 0„) with dimensions = A+ = 3 — A“, 
times ■ ■ ■ 77 ^". (For scalar helds heavier than 3H/2 the mapping is less clear since both 
asymptotic solutions, A^, decay the same way.) Now one may expect that the dominant 
squeezed limit behavior of correlation functions should be related to the OPE limit of the 
CFT operators, where two operators dual to the short modes are brought close to each other 
in position space, and their product is expanded as a sum over operators. We will show that 
this is not necessarily the case. 

As a concrete example consider (j){xi)(j){x 2 )(T{x 3 ) where 0 is a dimension zero held (an 
approximation to the almost massless huctuations of the inhaton). The coefficient of this 
term has the structure of (<h(a 3 i)<h(a 32 )S(a 33 )), where A$ = 3 and A^ = A/:. In the limit 
\xi 2 \ = \xi — X 2 \ ^ 0 , one expects 

01 ( 331 ) 02 ( 332 ) ~ ^cf2‘^a;f2‘^"'^^"'^''Q(33i), (3.10) 

Q 

where the constants cf^ are called OPE coefficients. Using {Q{xi)Q{x 2 )) oc 3312 ^'^'^ gives a 
term in log T of the form 

f Y\d^xi(j){xi)(j){x2)a{x3)x^^~^x^3^^ . (3.11) 

i 

Taking the Fourier transform, we obtain 

jd^kd^q(f){k)((){—k — q)a{q)f{k,q,ri), with f{k^q^rj)(xk^~^^q^^^~'^7j~^‘’. (3.12) 
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The 3-point correlation function can be computed via 

{<Pk,{v)<Pk2{v)(^k3{v)) = j T>0T>cT0(fci)0(fc2)a(fc3)|T[{(^},t]|l (3.13) 

For small non-Gaussianities, |Tp is dominated by the Gaussian part 


= exp 


W} 


(3.14) 


and the non-Gaussian correlation functions can be calculated perturbatively. The cubic term 
(3.12) gives the squeezed limit contribution 

{Mv)(l>-k-qiv)(^qiv))oPE = v)P}ik)Paiq, v) oc k^-^^P}{k), (3.15) 

where we have used the fact that 


Pa{q,v) tx (3.16) 

q-^ 

However, one can show that (3.15) is a subleading contribution in the squeezed limit. The 
leading contribution comes from correlating with Cglf]) a similar expression as (3.1) with the 
substitution 0. Up to an overall constant, this is fixed by the scale invariance 


(0fc0-fc-qO-g(h)) tx ?7 ^'^k P^{k)P^{q,ri). (3.17) 


(Note that super-horizon 0^ is time-independent in our approximation.) Hence 


i'Pk'P-k-q'^qiv)) OPE ^ ^ ^ ^ 

{(j)k(j)-k-qO'q{ri)) ^kJ 


(3.18) 


where we used the fact that A~ <3/2. Note that on the level of in-in diagrams the OPE 
result corresponds to having both power spectra on the short modes, while the diagrams for 
expression (3.17) have one of the power spectra on the long mode. These different diagrams 
are shown in Fig. 1. 

The dominant squeezed limit contribution (3.17) arises from the following cubic term in 
log T; 

J d^kd^q(j){k)(j){—k — q)(j{q)k~^’^ P^^{k)rj~^’^. (3.19) 

This is a local function of xis when transformed to the position space. Interpreted as 
a correlation function of dual operators, it is a contact term (<h(a;i)$(a; 2 )S(a; 3 ))^Qjj^g^^^ oc 
( 5 ^(a 3 i — 0 ^ 3 ). However, this does not make it any less relevant for cosmological correlation 
functions. We can translate to correlation functions of C, by using 


c ^ 


H(j) 

0s.r. ’ 


(3.20) 
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Figure 1: (a) The diagram that corresponds to the squeezed limit bispectrum calculated 
using OPE (see equation (3.15)). Both power spectra are on the short modes, (b) Diagrams 
for the leading part of the squeezed limit bispectrum given by (3.17). One of the power 
spectra is on the long mode. 


where 0s.r, is the background slow-roll rate of the change of 0. The held a would then 
contribute to the 3-point function if there is a mixing and, even if not, its exchange 

contributes to the squeezed 4-point function {\q\ = \ki + k 2 \ ki): 


(Cfci Cfc2 Cfcs Cfc4 ) (T 


H 


(0fci0fc2‘^q(h)) (0^3‘^’^4 *^9 (h)) 

Pa{q,v) 


(3.21) 


The resulting three and four point functions are non-analytic functions of q which lead to 
long-range correlations in position space. 

Our expansion (3.1) keeps only these local terms in logT, which as argued above play 
the dominant role in the squeezed limit. The corrections generically start at order estimated 
in (3.18). 


3.2 Entropic perturbations m/H 1 

The huctuations of a very light held a remain almost constant during inhation and the 
power in huctuations is nearly scale invariant. As a consequence they can ahect cosmo¬ 
logical observables at a later epoch of cosmic evolution, e.g. during reheating or when the 
a huctuations become time-dependent again (for instance when the expansion rate drops 
below their mass). In the latter case the mode functions are not given by (2.2) with A“ ~ 0 
anymore, however we will argue that the modes remain classical and to study the ahect of 
super-horizon perturbations of a on the observables one can still use an expansion in terms 
of local observables (3.1). 

Suppose at some time the solutions to the linear equation of motion for a become os¬ 
cillatory with a (time-dependent) physical frequency u. Now both mode functions have 
to be kept but the reason why the held can still be treated classically is that the huctua¬ 
tions have very large occupation number. The amplitude of the inhationary huctuations are 
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~ while the typical size of the vacuum fluctuations for a fluctuating held of 

frequency u is ~ l/a^u. The ratio of the two is a measure of the occupation number 

which is 


N 


JdL 

(^Dvac 


(k/a)^ 


(3.22) 


This is much larger than unity if the mode starts oscillating while it is still super-horizon, 
in which case k/a u Tfinf, and also when it enters the (sound) horizon long after the 
end of inhation, when a; = /c/a ~ TTreentry ^ Hini- This implies that a and a can both 
be known with great accuracy and the evolution of the held is determined in terms of its 
freeze-out value during inhation a. As long as the a mode is super-horizon this evolution 
is also local and therefore an expansion like (3.1) can be used to relate observables to the 
locally measurable operators made of a which can all be related to the local value of a and 
its derivatives. Moreover, since the interactions occur when all modes are already super¬ 
horizon they should be independent of k, therefore the squeezed limit behavior is the same 
as local non-Gaussianity. Note that, as before this local non-Gaussianity can be generated 
via a mixing between ( and or other composite operators made of a —energy density is 
perhaps the most natural candidate—while can in principle vanish. 


These models are known to violate standard single-held consistency conditions since 
unlike long wavelength ( perturbations those of a are not locally equivalent to a coordinate 
transformation. That is, while in the squeezed limit 


■Ccq , fcv) P {(ki ■ ■ ■ Cfcjv)) (3.23) 

where V is the dilatation operator, ■ ■ ■ , /cat) is completely model dependent. However, 

we can still hnd weaker conditions which are satished by models in which the violation of 
the single-held relations is triggered by mixing with other helds, thereby distinguishing them 
from models in which ( by itself evolves at super-horizon scales.^ 

The idea is to use the knowledge about the vertices with soft ( modes [like (3.23)], and 
the factorization to eliminate the ignorance about interactions of soft a modes. We hrst con¬ 
sider only a single entropic mode, which is easily generalizable to the case of several identical 
ones. Here, one can fully determine the degree of mixing between a and ( from combinations 
of observable squeezed-limit correlation functions. The mixing cannot exceed 1 by unitar- 
ity, otherwise the two-by-two correlation matrix will have negative eigenvalues. Moreover, 
diherent ways of determining it from diherent combinations of correlation functions must 
agree. This gives a consistency condition for these models. In the more general case, the 
unitarity requirement still puts a weaker constraint on the ratio of squeezed limit correlation 
functions. Most notably, on /atl and tatl, which characterize respectively the 3-point and 
the 4-point correlators when the sum of two momenta becomes soft. However, these weaker 
constraints on observables must be obeyed as a consequence of the axioms of probability 
[12], and therefore any consistent theory of fluctuations would automatically satisfy them. 

^The latter case can be characterized by the condition that the constant adiabatic mode of Weinberg [22] 
is a sub-dominant solution to the linearized field equations. 
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3.2.1 Consistency conditions on two-field scenarios 

Consider having only one entropy mode. In this case, we have 

iUM' = PcMUki) + (3.24) 


and 


iCkMM' =^dkiMk3)Pdq)+Uki)Uk3)Pa{q) 

+ ['Cc(^l)'Co-(^3) + ^a{kl)d{f^3)]P(a{q), 


(3.25) 


for the counter-collinear limit of the 4-point fnnction ki + k 2 = q ^ 0. Solving for r = 
P<:a/VPd^)Pdq): we obtain^ 

^ 2 ^ _ [iCkdk.Q' - n(g)ec(fci)][(a3a4C-.)^ - n(g)ec(fc3)] _ 

Pd^Mkdk^CM' - {Ckdk,Q'^df^s) - iCk^Ck^y ^df^i) + PdqM^iMf^dY 

(3.26) 

This must hold for all ki and k^ as long as they are much larger than q.^ Using the standard 
consistency condition (3.23), and assuming that r 1 simplifies the above equation to 

^2 ^ [iCkdk^Cd' - Pdq)PPdkimCk,CkX-d' - Pdq)PPdk3)] .3 27^ 

PdQ){Ckdk,CM'-P!iq)PPdki)PPdk3) 


The above equation tells us that, even though the a field itself is unobservable, the 
mixing parameter r can still be measured just from the observable correlation functions 
of (. Moreover, r can also be solved in terms of the counter-collinear limit of the 5-point 
function and the squeezed limit of 3- and 4-point functions, and so on. Then, given that 
the two expressions for r must agree, one can write down an equality that involves just the 
correlators of (. These relations are new consistency conditions for the two-field inflation 
that hold irrespectively of the details of a interactions. 


3.2.2 Unitarity constraints on multi-field models 

In the more general case where a set of local operators {Oa} contribute to the squeezed limits, 
the constraints are weaker, since in general one cannot solve for the individual mixings r^. 
However, there are inequalities that must be satisfied for the theory to be unitary.® We 
derive the simplest one: The squeezed limit of the 3-point function and the counter-collinear 
limit of 4-point function now read 

iUAY = J^PcMUkd + PdqM^i), (3.28) 


^Notice that in this equation we treat r non-perturbatively. 

^One may wonder if the denominator can be equal to zero for some momentum configurations. It is equal 
to Pa(<z)^’c('?)^o-(^i)^<T(^3) therefore it can vanish only when the numerator vanishes as well. Note that 
in this case r cannot be determined using (3.26). 

®For a discussion of unitarity constraints on single field inflation see [23]. 
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and 


{Ck,Ck,Ck,Ck,y =ec(fci)ec(fc3)Pc(g) + J2Uki)Uk3)Pab{q) 

a,b 

(3.29) 

+ + ^a{kl)^({k3)]P(^a{q), 

a 

where Pat = {OaObj- Choosing an orthogonal basis for the operators such that Pat oc 5ab, 

and assuming for simplicity that ki ~ k^, we arrive at the following constraint (using the 
Cauchy-Schwartz theorem) 

[(a.a.c.)'-n(?)ec(fci)]^< (E^')^c(?)[(a,Cfea3a4)'-2(a.a.c.)'ec(fcl)+n(?)ec(fcl)ec(fcl)]• 

(3.30) 

In the extremely squeezed limit where only local shape survives, this simplihes to 

(I/™)' < (Z-'S - ((I- (I/“)')] • p.31) 

where f^i is the observable part of the local non-Gaussianity: the difference between to¬ 
tal Jnl and the single-held part. A stronger bound can be obtained by subtracting the 
contribution of tensor mode exchange from the 4-point function since = 0. 

On the other hand, 1 — is proportional to the determinant of the total 2-point 

correlation matrix detPi^ip/, where ip indicates all local helds and operators including (. The 
positivity of this determinant together with (3.31) then implies 

<tnl- (3.32) 

* * * 

Let us close this section by two general remarks. First, the underlying assumption for 
the expansion (3.1) was that the short modes are in vacuum until the long mode has crossed 
the horizon and become classical. This assumption can be broken for instance if the short 
wavelength modes get excited inside the horizon. An explicit example of this sort is axion 
models of inflation [24, 25] where there is another physical frequency ca > P in the problem 
at which the modes get excited. This case has been studied in appendix A.4. As seen after 
sufficient squeezing (here q/k H/oj) our formalism will be applicable again. In particular, 
in the absence of other degrees of freedom the single-held consistency conditions hold for 
sufficiently squeezed correlators [26]. 

The second remark is regarding the resemblance between the expansion (3.1) for classical 
long wavelength modes and the bias expansion for the distribution of halos and galaxies. In 
the latter case the dynamics of the short scale modes is inhuenced by the local observables 
of the long wavelength perturbations of gravitational potential (made of didj(j)) along the 
trajectory of the short modes. However unlike a super-horizon a held whose time-evolution 
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is local and fully determined in terms of the field a and its spatial derivatives at a single 
time-slice, the sub-horizon density fluctuations have long range gravitational interactions. 
Therefore the bias expansion won’t be local if written just in terms of the initial held and 
its spatial derivatives (see e.g. [27]). 


4 Heavy fields m/H > 3/2 

In this case the two solutions of the free-held equation decay with the same rate and oscillate, 
scaling as ^ — 9/4. Hence, the helds are not classical and we 

cannot expand as in (3.1). However, there is still a sense in which the correlation functions 
arising from the exchange of heavy helds factorize in the squeezed limit. This is because in 
this limit the dominant contribution arises from the production of on-shell heavy particles, 
whose amplitude is suppressed by a Boltzmann factor exp(— tt/x) for each pair of particles. In 
the weakly coupled case the other contributions in which no on-shell particle is created are 
captured by a set of ehective derivative self-couplings of ( arising from integrating out the 
heavy particle. These self couplings are not suppressed by any Boltzmann factor, but they 
lead to equilateral type non-Gaussianities. For instance, the squeezed limit bispectrum scales 
as which decays by a factor of {q/kY^'^ faster compared to the former contribution from 
production of on-shell particles which goes as {q/kY^"^. For higher statistics these on-shell 
particles give a sub-dominant piece but they still have a characteristically different functional 
form. This problem has been thoroughly studied in [13, 14] and we will just touch upon the 
salient features of the squeezed limit correlators.^ 

Instead of correlation functions in a classical background, here we should talk about the 
wavefunction of universe to keep track of the phase of the a held. Since the concept of particle 
is well-dehned for late-time heavy helds, it is more natural to expand the wavefunction in 
terms of eigenstates of particle number operator as in the scattering theory. Let us choose an 
intermediate time fj satisfying k~^ —fj q~^, so that the super-horizon held aq is well- 
dehned but the hard k modes are in vacuum. The wavefunction at ? 7 o, when all correlations 
are formed, can be written as 


l'l'(%)) = E /n^ /DC u(, |{C}.iV^> l-KO.iVaiij], (4,1) 

N n=l ^ 

where U{r]o, fj) is the time-evolution operator and the sum is only over iV-particle states of a 
with super-horizon wavelengths. Other possible degrees of freedom have been omitted. The 
wavefunction |4'(i7)) depends on the sub-horizon dynamics, and if a exists there as a weakly 
coupled degree of freedom can be perturbatively calculated. However, on general grounds 
we expect that 

(4,2) 

since we assume that the theory is in the adiabatic vacuum deep inside the horizon, and 
hence particle production is a result of evolution in the expanding background. Secondly, 

^For a study of the impact of heavy fields on correlators involving tensor modes see [28]. 
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all modes that are still deeply sub-horizon must be in vacuum; 

'h[{C}, ^0-; i?] oc (5(Cfc), for all —/cr^ 3> 1. (4.3) 

(The ie prescription evolves such a wavefunction into the true vacuum.) On the other hand, 
the soft (g modes are constant up to corrections of order (gr))^, therefore 

(C,|t/(%, ijjlC,) « ^IC, - C, + 0(ef)\. (4.4) 

The correlation functions are obtained by taking averages using the wavefunction (4.1). 
Keeping the minimum number of Boltzmann factors, we get: 

(CqCfciCfea),, ^ (CfeiCfcallc^q) J dCg Cg ¥[Cg,Oa;fi] T[Cq,la_q;r)] + C.C. (4.5) 

The second factor on the r.h.s. depends on the sub-horizon mixing of (g and ag. We expect 
from the perturbative calculation that this correlation should form at —qri ~ ^ 1 and 

to introduce another factor of exp(—/i7r/2). The first factor on the r.h.s. corresponds to a 
massive particle decaying into a pair of (k modes when their total energy red-shifts to values 
around the mass, i.e. r] ~ —^/k. The super-horizon wavefunction Ug oc c.c. 

then leads to an oscillating piece cos(/ilog/c -1- (pg) and a dilution factor (g/A;)^/^. The phase 
4>g is fixed if we require that a should be in the adiabatic vacuum at very early times. 
As emphasized in [14] the oscillations result from interference: the expectation value is 
calculated between the Gaussian part of the wavefunction, and the part with production of 
a pair of entangled a particles, one of them oscillating into a (g mode, and the other decaying 
at a much later time into a pair of short wavelength modes. Fixing the other factors by 
approximate scale invariance, we obtain 

(C<?CfeiCfc 2 )a 0^ Pd(l)Pdk) (I) ' cos(/ilog| - 0o)- (4.6) 

This is enhanced in the sqneezed limit by {klqY^^ compared to equilateral type shapes. 
Moreover, the oscillating component makes this shape distinct from other contributions. We 
could also consider the a-exchange contribution to the 4-point correlator: 

dkiCk2Ck^Ckdl a OC ^ ^ (Cfci Cfc 2 I | l(J_q,) -\- C.C. (4-7) 

Again this is an interference between the Gaussian piece, and a piece which describes the fol¬ 
lowing process: (i) An entangled pair of massive particles is prodnced in the time-dependent 
backgronnd. (ii) The particles move ont of cansal contact, (hi) They subseqnently decay 
into several C, helds, when the total energy of the red-shifting modes coincides with m. 

Finally, if the exchanged held has spin s > 0, the sqneezed limit resnlt would be propor¬ 
tional to ■ ■ ki^. Only the longitudinal polarization of a massive tensor held can 

mix with a scalar e.g. via C,di^ ■ ■ ■ Therefore, 

(CqCfciCfc 2 )^ oc PYq)Pdk) (I) ' cos(/ilog| - (j)o)^s{q ■ k), (4.8) 

where is the Legendre polynomial. 


16 



4.1 Quantum versus classical correlation 

We have seen that the squeezed limit correlation functions have an expansion in terms of a 
discrete set of power-laws 

(c,&.a,) ~ n(9)-Pc('=)E“^ ilT • 

^ A (4-9) 

(CfciCfc 2 Cfc 3 Cfe 4 ) ~ Pd^i)Pd^3) ^ j counter-collinear , 

A 

with A integers > 2 for equilateral type non-Gaussianities, real numbers corresponding to 
the light {m/H < 3/2) degrees of freedom and their derivatives and products, or pairs of 
complex numbers corresponding to heavy degrees of freedom. This expansion is unambiguous 
and non-trivial and seems to be dictated just by the isometries of a quasi-de Sitter space- 
time. As emphasized in [14] the case of heavy held exchange is fundamentally different as 
it corresponds to a quantum mechanical interference effect.® In fact one can see that the 
fundamental difference between the case of heavy intermediate helds and the other cases 
is the difference between quantum versus classical correlations. The outcome of two short 
distance measurements in two far separated laboratories are correlated with each other in 
all cases, but in one case they are deterministic while in the other case the outcome of one 
measurement inhuences the result of the other. 

This difference has a distinct imprint in the squeezed limit correlators: When the long 
wavelength mode becomes classical, the a held works as a hidden variable and the expansion 
(3.1) is possible. As a result the dependence on the momenta on two sides of the intermediate 
long wavelength mode factorizes (see for instance (3.3) and (3.4)), or more generally it can be 
written as a sum of factorized terms weighted by the correlations among various derivatives 
and products of a. Different short distance measurements have no inhuence on one another; 
they are correlated because they are both inhuenced by a single long mode. On the other 
hand, the heavy held contribution is not factorizable because of the dependence of the 
squeezed 3-point and 4-point functions, respectively, on \og{q/k) and log(g^/A;i 2 /c 34 ) (where 
kij = \ki- kj\) [14]. 

5 Prospects in large-scale structure surveys 

In this section we will discuss the observability of various squeezed limit behaviors, when 
contrasted with the contribution from the gravitational interactions during the structure 
formation (see [26] for a related analysis in the context of single-held inhation). Consider a 
non-Gaussian 3-point correlator whose squeezed limit behavior is 

(C,a4') = fNLP(k)P{il) + k + k’). (5.1) 

®It was also argued that they are distinguishable due to the non-locality of the real space correlator. 
However, almost all terms in the above expansion lead to large distance (non-local) correlations when Fourier 
transformed with respect to q (see the discussion of section 3.1.1). 
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The signal to noise ratio in a 3d survey of volume V can be estimated as 


(S/Nf 


(27r)9 


d^kid^k' d?q 


{CqCkCk') {CqCkCk') 
{CqCkCk'CqCkCk') 


(5.2) 


Using (5.1), approximating the denominator by the Gaussian contribution, and replacing 
(27r)^(5(0) —)■ V (appropriate for hnite volume surveys), we get 




(5.3) 


where we are calculating the signal in the squeezed limit configurations q k. The k 
integral is dominated by the highest values of k. This is generally the case since the number 
of modes grows as k^. For biased tracers, k really corresponds to the modes that have become 
nonlinear. The dependence on k would then be different. However, we are interested in the 
dependence on q as it is made small so that (q mode is linear and P{q) oc 1/q^. Limiting 
the inner integral to g < q^ax we get a contribution which is proportional to 


(S//V) K qt 


(5.4) 


That is, the error with which this type of non-Gaussianity can be measured by considering 
modes longer than scales as This suggests that unless A = 0, there is more signal 
in larger values of g. 

However, there is a guaranteed level of non-Gaussianity due to the subsequent grav¬ 
itational interaction of the modes that also increases as g becomes bigger and competes 
with primordial signal. For very large values of the hard momenta k, this is not accurately 
calculable. However, in the squeezed limit its g dependence must scale as 

iCqCkCk'),, OC 0(g')P(g), (5.5) 

since the short scale dynamics can only depend on the local observables made of the long 
mode. As k is kept fixed around the largest available values and g is made smaller, this 
gravitational contribution decays faster than (5.1) for all A < 2. Therefore, the relevant 
question is how the signal in modes with g < gmax, he. equation (5.4), compares to the 
ratio of (5.1) to (5.5), which scales as q ^~^ . One hnds that for A < 1 the primordial shape 
is better distinguishable in the squeezed limit, while for higher values of A analyzing the 
equilateral conhgurations signihcantly improves the chance of detectability. However, unlike 
the squeezed limit the equilateral primordial shapes are not expected to be very distinct 
from gravitational contributions, hence the theoretical uncertainty in determining the latter 
close to the nonlinear scale poses a challenge to increasing k. 


6 Conclusions 

We argued that several squeezed limit properties of non-single-held inhationary models hold 
more generically, and can be studied in a unihed fashion. The helds in a quasi-de Sitter 
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space-time become free at super-horizon scales and follow power-law time evolution. The 
light fields become classical and the leading squeezed limit behaviors can be obtained from 
an effective parameterization of short-distance correlations in terms of a local expansion in 
derivatives and powers of the background long wavelength mode. We showed how the known 
results, such as the squeezed limit behavior in quasi-single-held models, and the unitarity 
constrains on mnlti-held models, can be derived in a more model independent way from 
this formalism. Heavy helds, on the other hand, have complex wavefnnction at late times 
and are natnrally described in terms of many particle states. The sqneezed limit correlation 
fnnctions due to the production of these heavy particles by the time-dependent backgronnd 
is dominated by interference effects and hence they oscillate as a function of the ratio of 
kshovt/k\ong- All these features have not much to do with the details of sub-horizon dynamics, 
but they essentially depend on quasi-de Sitter evolution and the existence of (not too heavy) 
snper-horizon degrees of freedom other than metric and inhaton huctnations. The form can 
be used as a general test of a quasi-de Sitter cosmic evolution. Finally, we discussed the 
observability of various squeezed limit scalings of the bispectrum. 
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A Purturbative examples of factorization 

A.l Fundamental fields 

Suppose C is coupled to a light field a via 

H=Hi + H2 = CV + Cct. (A.l) 

We are interested in the leading contribution of a u exchange to the squeezed limit of the 
3-point function (CqCfciCfca)^- 

(CqCfciCfc2)^2 =-2Re/Cq(r?o)Cfci(ho)Cfc2(ho) / ^^i(hi)^2(h2)\ 

\ J —OQ Vl J — OO V 2 / 

(A.2) 

+ / / ^Wi(r;i) Cg('to)Ci,i(fto)Cte('to) [ ^^ 2 ( 1 ) 2 )) + (-» W2(r)2) , 

\J—OO Vl J—OO V2 / 

where we introduced an infrared regnlator tjq since the held a is generically time-dependent 
at snper-horizon scales. We will often drop the time argnment of ( if it is r]Q. The main 
contribntion to the squeezed limit is when the mixing 712 between the soft ( and a helds 
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occurs at a time around rjg ~ —^/q much earlier than the "Hi interaction at r]k ~ —1/k. Sub¬ 
stituting (A.l) in (A.2), expanding in terms of mode functions and using this approximation, 
we obtain 


-Cq(ho)Cfci(^o)Cfc 2 (ho) [ ^ [ ^CkAVi)CkiVl)^q{Vl)(^*q{V2)Cq{V2) + C.C. 

J-oo 71 7-00 72 

+ CqiVoXl.iVoXkiVo) [ ^CkAVl)Ck2{Vl)(^q{Vl) [ ^(^q{V2)Cq{V2) + C.C. 

J —oo m j- oo Vi 

We next argue that this expression factorizes into 


{CqCkiCk2)a- ~ 


(Cq^-q(ho)) (CfclCfc2^q(h0)) 

Pa{q,Vo) 


oc P{q)P{k) I 


qy 


k) 


(A.3) 


(A.4) 


where A^. = |~Y|~^is the scaling of the growing mode and Pa{q,Vo) = Wq{Vo)\‘^- The 
two mixed correlators in the numerator are given by 

po 

{Cqa-q{rio)) =-iCq{rio)a-q{rio) -^(^-q{v)Cq{v) + c-c. (A.5) 

J — oo V 

and 

dn 

(CfciCfe2C7q(ho)) = -<fci(ho)Cfc2(ho)f^<7(ho) / ^Cfei(hl)Cfc2(hlX(hl) +C.C. (A.6) 

7-00 Vl 

The relation (A.4) is then verified by noticing that: (a) The integral in (A.5) is independent 
of ? 7 o as long as —gr/o ^ 1- Therefore, since the rji integral on the hrst line of (A.3) is 
dominated by —rji ~i//c<^l/g^we can replace the upper-bound of the 172 integral with 770 . 
(b) At late times the time-dependence of the mode functions factorizes and we can replace: 

c7g(r/2) ^a*(r/2)^^. (A.7) 

^q\VOj 


A. 1.1 Light fields 


For helds of dimension A^- < 1, the integral in (A.5) is naively dominated by late times, and 
the approximation (a) seems to be invalid. However, this late-time dominance is spurious 
since it is pure imaginary and cancels in the final result. This independence from the late¬ 
time cutoff 770 (apart from the trivial scaling in cTg(77o)) becomes manifest by reorganizing 
the perturbation theory as time-evolution problem [29]. In this language there are two 
contributions to (CqCT_q( 770 )): 


{Cq(y-q{Vo)) = (^-qivo) 


^^0 rj^ . 

-^Gt{VO,V)(^-qiv) + CqiVo) 

-00 V 


dv^R 


.(ho,h)C(h), (A.8) 


where 

Gf^iVop) = ^(/q(ho)/q(h) - c.c.) 


(A.9) 
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is the retarded propagator of the indicated held. At late times it decays as while (qiv) oc 
rf and (JqiTj) oc . Hence, for all A^- > 0 each term in (A.8) is manifestly IR-safe, i.e. it 
does not depend on the upper limit of the integral as long as —qrjQ 1, but only the sum 
of the two terms is real. The case of an unprotected massless held A^- = 0 is special. Here, 
the dependence on ? 7 o is real and corresponds to the possibility of the generation of a mass 
for a via its mixing with C,. 


In this formalism, the leading contribution of a exchange to the squeezed limit comes 
from two diagrams 


r>J?o 


(CqCfciCfea) = / ^Gl{r] 0 ,r] 2 )a-q{ 7 ] 2 ) [ 

J-oo 72 


'-oo Vi 


r‘??o 


dr] 


I'm 


dV: 


(A. 10) 


J —OO 'll J —oo 72 


where dots represent the part of the diagrams corresponding to the interaction among the 
high momentum modes at rji ~ —1/k. The approximation (b) can now be rephrased as 


G^^{r]i,r]2) G^^{r]o,ri2) 

and leads to the same factorization formula (A.4). 


for 


qVi < 


(A.ll) 


A.2 Local operator exchange 


Suppose C is coupled to a light held a of mass m < 3H/2 via 


n = ni + H2 = + Ca^. 


(A.12) 


We are interested in the leading contribution of a a loop to the squeezed limit of the 3-point 
function (CqCfciCfca)^- 


(CqCfciCfca)^^ =-2Re/Cq(r?o)Cfci(ho)Cfc2(ho) / ^'Hi{r]i)n2{v2)) 

\ J — OO Vl J —oo ^2 / 

(A.13) 

+Re/ f %-'Hi(ih)Co('?i))C/.,(%)Cfe(%) / ^W2fe)')+'Hi(>?i) W2fe) , 

\J — oo Vl J —oo V 2 / 

where we introduced an infrared regulator 770 since the held a is generically time-dependent at 
super-horizon scales. We will often drop the time argument of C if it is 770 - The contribution 
to the squeezed limit is insignihcant unless the momenta of the virtual pair and ^2 of 
the same order as the soft momentum q, when an 7/2 interaction among the three soft modes 
occur around ~ —i/q much earlier than the "Hi interaction at 77 ^, ~ —1/k. Substituting 
(A.12) in (A.13), expanding in terms of mode functions and using this approximation, we 
obtain 


cm 


-2Re 


d77i n° dr ]2 


Roo Vi J_oo V 2 L 


CqiVoXki {Vo)Ck2 (ho)Cfci (hl)Cfc2 iVl)Cq (^ 2 ) 


(A.14) 


-Cq(h0)Cfci(h0)Cfc2(h0)Cfci(7/l)Cfc2(7/l)Cq(h2) {(^l{Vl)(^-q{V2)) , 
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where 

Using the late time behavior: 


(Tq, {ril)(^q^ (h2)- 


'<5fl+<?2=<5f 


<32' 


(A.15) 




_ 3 /9 m2 

" “ 2 “ Vi “ hP 


(A.16) 


in (A.15) and substituting in (A. 14), it is easy to see that it factorizes: 

{Cqa_q^{riQ)a_q^{rio)) {Ck^Ck2(^q,{Vo)(^q2iVo)) 


{CqCkiCk2) rr — 


'<3i+<32=<3 


PAqi,vii)PA<i2,vii) 


(A.17) 


However, one can show that the above squeezed limit 3-point function fully factorizes. Note 
hrst that because of the late time behavior (A.16) 


CkiCk2(^q,{Vo)(^q2{Vo)) = f{ki,r]o)P^{qi,rio)P^{q2,rio) , 


which implies 


{CqCk^Ck2)^2 ^ f{kl,V0) {Cq(^-q{V0)) ■ 
Next, integrate (A. 18) over and q 2 , with qi + q 2 = Q, to obtain 

{CkiCk2(^l{Vo)) = f{ki,rio)P^2{q,rio). 
Combining (A. 19) and (A.20) we hnally obtain 

■..a,'..I,,. 


(A. 18) 
(A. 19) 

(A.20) 

(A.21) 


Let us make a few remarks: 


• The above derivation holds only if A^- < 1 . Otherwise, the r]i integral in (A.13) is 
dominated by r]i ~ 772 and not by ? 7 i ~ — 1/g: the UV divergence in the a loop wins 
over slow variation of Cij(hi) ^t super-horizon scales. Nevertheless, by renormalizing 
the composite operator one expects to obtain the same factorized contribution. 

• Although the above formula was derived for the operator it is expected to hold 
more generally for any intermediate operator O 




{Oq{r]o)0_q{r]o)) 
and the last scaling to be valid even for large mixing. 
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A.3 Soft internal lines 


The approximation (b) of section A.l can be used to also show that the contribution of soft 
internal lines (or operators in general) to the correlation functions factorizes. For instance, 
a 4-point function in the counter-collinear limit ki + k 2 = q ^ 0 can be written as 


{CkiCk2Ck3Ck4) o ~ 


(CfclCfc 2 ^< 5 f(ho)) (Cfc3Cfe4f^-(j(ho)) 
{Oq{rio)0_g{rio)) 


(A.23) 


The intermediate operator can now be a soft (q mode. 


A.4 Resonant non-Gaussianity 

The above predictions always require sufficient amount of squeezing q/k 1. This is needed 
to ensure that the short modes are in vacuum before the long mode crosses the horizon. An 
illustrative example is resonant non-Gaussianity where the factorization on soft internal lines, 
e.g. (A.23), can fail due to insufficient squeezing. Here the modes get excited at physical 
frequencies of order the resonant frequency u ^ H and the interactions start to be important 
when the short modes cross this threshold. If q/k > H/u the soft intermediate mode is not 
yet super-horizon when the interactions of the short modes become important. Thus, the 
approximation (b) cannot be made. More explicitly, the leading contribution to a 4-point 
function in the counter-collinear limit contains a term 

-Cg(ho)Cfci(ho)Cfc 2 (ho) [ 

J-oo vt J-oo vt (A.24) 

Cfci(hl)Cfc2(^l)Cq(^l)Cfe3(h2)Cfc4(^2)Cq(h2) + {^1,^2} ^ {^3,^4} + C.C. 

where V^'^\r]) oc cos ut is the leading cubic vertex. The integrals are dominated when 
—krj u/H around the saddle point. Suppose q k but q/k ^ H/u then we can neglect 
the factor exp(—ig?]) in Cg(h) replace 

Cqiv) = Cq(ho)(l + ^ iqvCqiVo), (A.25) 

which gives an opposite contribution compared to Cq{v)Cq{Vo)/Cq{Vo) fhe same approxi¬ 
mation. Hence the approximation (b) does not hold in this regime. 
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